Potential-density pairs for a family of finite disks 
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ABSTRACT 

Exact analytical solutions are given for the three finite disks with surface density 
S„ = o"o(l — Ft? / a 2 ) n ~ l l 2 with n = 0,1,2. Closed-form solutions in cylindrical co- 
ordinates are given using only elementary functions for the potential and for the gravi- 
tational field of each of the disks. 

The n=0 disk is the flattened homeoid for which T,h om 



cro/y 1 — R 2 /a 2 . Improved 
results are presented for this disk. The n=l disk is the Maclaurin disk for which 
^Mac = o~q^/1 — R? /a 2 . The Maclaurin disk is a limiting case of the Maclaurin spheroid. 
The potential of the Maclaurin disk is found here by integrating the potential of the 
n=0 disk over a, exploiting the linearity of Poisson's equation. The n=2 disk has the 



surface density = <?"o (l 
potential of the n=l disk. 



R 2 /a 



2\3/2 



The potential is found by integrating the 



Subject headings: gravitation- methods: analytical-celestial mechanics- galaxies: kine- 
matics and dynamics - galaxies: individual(NGC 891) 



1. Introduction 



X 



Observations of edge-on galaxies now provide detailed structural and kinematic 3-D informa- 
tion. However attempts to use this data to generate mass distribution have not been completely 
successful. One problem is that it is difficult to calculate the potential and force vectors for finite 
disks. There is a need for fully solved finite disks which can be used for theoretical studies, for 
benchmarking computer programs, and as basis functions to directly model 3-D observations. 

Mass modeling commonl y assumes that the disk mass in contained in an infinite disk. These 
include the exponential disk (Freeman 1970), the Mestel disk (Mcstcl 1963), the Kuzmin-Toomre 



disk (Toomre 



1963 



Rvbicki disk (|Evans k Collettlll993h 



Bin nev fc Tremaind 12008 : Evans fc de Zeeuwl 1 19921 : IConwavl l2000f ). and the 



Only a few finite disks h ave been solved analytically for all (R,z). lLass fc Blitzed ()1983l ) and 
Vokrouhlicky fc Karad (|199 8) give a solution for a for all (R, z) for a finite thin disk with constant 
density. The gravitational attraction approaches infinity at the edge of this disk. The gravitational 
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attraction of the finite Mestel disk (IMestellll963 ; iLvnden-Bell & PineaultJll978l ; lHunter et al.lll984l ) 
and the truncated exponential disk (jCasertand Il983l ) are finite a t the edge of the di sk bu t these 
disks have not been solved in closed form for points off the disk. iHure et al.l (|2008l ); iHurel (j_2005) 
describe a method to approximate the potential of a power law disk for points both on and off the 
disk. 

The family of finite disks with surface density T, n (R; a) = (l — R 2 /a 2 ) n are related to the 
Maclaurin sphe roid which has been stud ied s ince the time o f New t on. This family h a s rece ntly 
been studied by Gonzalez h Reina ( 2006 ) and Pedraza et al.1 ( 2008 ). Gonzalez k, Reina ( 2006 ) use 
the method of iHunterl (|1963l ) to obtain the general solution as a sum of Legendre polynomials 
in elliptical coordinates and give evaluated expressions for the potentials for the disks 1, 2, and 
3. Here we derive complete closed form solutions in cylindrical co-ordinates for the potential and 
the gravitational fields of the n=0, 1, and 2 disks. We simplify the integration by moving to the 
imaginary domain in way which is similar to the com plex-shift method introduced by Appell. See 
Ciotti &: Marinaccil (J2008J) ; ICiotti Sz Giampieril (|2007l ) and references therein 



The gravitational potential follows Poisson's equation V 2 $ = 4nGp. Poissons's equation is 
linear so that solutions may be added. That is, if V 2 ^! = AirGpi and V 2< I>2 = ^7rGp2 then 
V 2 (<I>i + $2) = 4:irG(pi + p ?). Simil a rly so lutio ns may be dif f erentia ted with respect to a parameter 



to obtain a new solution. 



Toomre! (|l963h and iLvnden-BelJ dl989h exp loited this linearity to find 
fam il y of v elocity-density pairs by differ entiating thelKuzminl (1956 ) model . In the same way, 



Satohl ()1980l ) derived new models from the Miyamoto k, Nagail (119751 ) set of potential-density pairs 
by differentiating with respect to a parameter. Here we use a similar approach by integrating 
a known solution with respect to a parameter. If the integral is tractable the result is a new 
potential-density pair. 

Maple 11 was used in this work. Maple was invaluable in simplifying the unwieldy intermediate 
results but needed a good deal of coaxing to handle the messier equations. The Maple results were 
cross checked in several ways. 

A common notation is used throughout. Cylindrical co-ordinates (R, z) are used; £ is the disk 
surface density; <jq is the surface density at R = 0; and a is the disk radius. $ is the potential 
where $ is always negative and $(00) = 0; Fr and F z are the gravitational field vectors with the 
standard sign convention, ie: F = — V<£ . 

The principal values of elementary functions are used so that, for instance, for z = x + I y = 
rcos(6*) + I r sm(9); yfz = -y/r cos(#/2) + I sin(#/2), valid for — ir < 9 < tt. In this way the 
function y/z is unambiguous with continuous derivatives except on the negative real axis where it 
is discontinuous. 



- 3 - 



The n=0 Disk: The Flattened Homeoidal Shell 



A homeoid is a shell of uniform density which is bounded by similar spheroids. Newton was 
the first to prove tha t the net force is (i.e., the potential is constant) within these shells. See 
Chandrasekharl (jl987l ) for historical background. The infinitely thin homeoidal shell is a differential 
element of a spheroid. The homeoidal disk is the limiting case for which the minor axis approaches 
0. The collapsed density of the thin homeoid is: 



J hom 



(R;a) 







l-R 2 /a 2 forR<a 
for R > a 



Lvnden-Belll (|1989l ) gives a formula for the external potential of the thin homeoidal shell. 
Taking the limit c = (1 - e 2 ) l / 2 a -> gives a solution for the disk which is valid at all R and z. 



Cuddefordl (|1993l ) gives an expression for the potential of this disk which is simpler than 



previous solutions: 

&hom(R, z; a) = — 27rao"oGsin _1 



2a 



^z 2 + (R + a) 2 + ^z 2 + (R - a) 2 
Equation [2] can be further simplified. First make the trivial transformation: 

yjz 2 + {R + a) 2 - yjz 2 + (R- a) 2 



(2) 



®hom(R, z;a) = -27rao- Gsin _1 

Now use the identity I All to obtain 

® hom (R,z;a) = -iraaoG 



2R 



(3) 



sin 



-i 



a — I z 
R 



+ sin 



-i 



a + I z 
R 



(4) 



Equation[4]can be integrated over a whereas Equation [3] yields an impossible integral. Equation 
Hlmust be real valued on physical grounds. This is easy to prove by noting that sin _1 (x) is an odd 
function of x and so the odd powers of I z in the series expansion of equation H] cancel, leaving a 
real valued result. 

The gravitational field for the collapsed homeoid is obtained from the potential. Equation [4] 
yields particularly simple expressions for the field vectors F^om and F z ^ om : 



FR,hom(R, z; a) 



iraaoG 
R~ 



a — I z a + I z 



Fz,hom(R,z;a) = -Traa G 



{a-Iz) 2 ^R 2 + (a — I z) 2 
I I 



^R 2 -(a -I z) 2 ^R 2 - (a + I z) 2 



(5) 
(6) 
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These force vectors can be expressed as entirely real functions using identities IA8I and [ 



Where 



F RMm (R, z; a) = -V27rao- ^ 

Rjih 

77i / 73 \ ^ sgn(z) v / /i/2 + /3 
F Z}hom (R,z;a) = -V2iraa ( J — 

/1/2 



(7) 

(8) 



/i = ^/z 2 + (R + a) 2 

h = ^z 2 + (R-a) 2 (9) 
h = a 2 -R 2 - z 2 



The n=l Disk: The Maclaurin disk 



Beginning in the early 18'th century Co lin Maclaurin, along w ith James Ivory and many others, 
studied the proper ties of elliptical bodi e s. IChandrasekharl (119871') includes a v ery good historical 



(119681 ): iKalnaisI (I197ll . Il972 l. 



summ ary. See a lso binnev Tremainel (l2008h : bertinl ^OOol ): ISchmidtl |l95^); iMihalas Routlvi 



The homogeneous oblate spheroid is the simplest case of a spinning body for which the gravi- 
tatio nal attraction balances the centrifugal force. The Maclaurin disk, also known as the Kalnajs 
disk (|Kalnaislll972l ). is a limiting case for which minor axis is 0. The Maclaurin disk is defined by 
the surface density: 

for R < a 
for R > a 



-■Mac 



{R;a) 




R 2 /a 2 



(10) 



There are a few solutions for the potential of the Maclaurin disk in the literature. IMihalas Routlv 



(jl968l ) gives an e xpres sion for the potential of an oblate homogeneous spheroid based on the deriva- 



tion in ISchmidtl (|1956l ). The potential of t he Macl a urin d isk can be found by letting the eccentricity 
e — > 1 while holding the mass constant. iHunterl (|1963l ) give s the solution for the Maclaurin disk 
as a s e ries of Legendre po l ynom ials in elliptical coordinates. iNeugebauer fc Meinell (|1995l ): iMeinel 
(|200ll ) ; iGonzalez &: Reinal (|2006l ) give a closed form solution for the potential of the Maclaurin disk 
in elliptic coordinates. 

The starting point here is the potential-density pair for the n=0 disk, the flattened homeoid 
for which E/,, om (i?; a) = cro/i/l ~~ R 2 jot 2 . The surface mass density of the Maclaurin disk is found 
from the transformation: 



^Mac{R;a) = - I T, hom (R;a) da = — [ - °"° = 
a Jo aJoy/l- R 2 /a 2 



da = (Tq yl — R 2 /a 2 



(11) 
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The corresponding potential is: 



1 f a r -7R7 G 

®Mac{R, z;a) = - / $hom{R; a) da = 

a Jo a 



a — I z 
R 



a + I z 
R 



n Mil ' | — | — si u ' | — | da 

(12) 

Where the expression for $fe om is given by equation H] above. Use integral 2.813 and 2.833 from 
Gradshtevn k Rvzhikl (|l99J) to obtain: 



$Mac(R,z;a) 



7ro"oG 
4a 



(2af - Bf + 2z z ){ sin 



1 fa + I z 
B 



+ sin 



a — I z 
B 



+ a y/R 2 - (a - I z) 2 + yjR 2 - (a + I z) 



(13) 



3 z lV^ 2 -(a-I^) 2 -Iv / ^ 2 -(« + 1 ^) 2 



Equation [13] can be converted to an entirely real expression by using the identities IA11 IA6l and 



$Mac(R,Z;a) 



7T(ToG 

4a 



2(2a 2 - R 2 + 2z 2 )sin _ 



i / A — h 
2R 



+ - / 3 - 3^2 | z | y/Kh + h 



(14) 



where /i, /2, /3 are given by equation [9] above. 

The gravitational field of the Maclaurin disk is found from the potential using $ as given 
by equation [T3l The resulting expressions were expressed as entirely real functions by using the 
identities EH \M[ and \M\ : 



FR,Mac(R,z;a) 



7TCJoG 

2Ra 



,, ,,2 ■ -l ( h~ h 



AT / 2 r>2 i 2\ V /l/2 — /3 

+ V2a(a -it + z ) — 

hi 2 

AT 1 1 /- 2 , D 2 , 2n v 7 /1/2 + /: 
- V2 I z I (a +-R + z ) — 

Jl/2 



(15) 



.Mac 



(i?, z; a) 



7T(7oG 



a 



2V2. 



az- 



2zsin 



shh — h 



-1 ( fx — /2 



2it 



/1/2 



AT m / 2 d2 2^/1/2 + /3 

+ V2sgn(z) (a - it -z) — 

J1J2 



(16) 
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where /i, /2, f3 are given by equation [9] above. 

The potential on the z axis and on the z = plane can be found by by taking limits of equation 

7TO"oG 



a 



a 



+ z 2 ) sin -1 



a 



\J a 2 + z 2 



a z 



(17) 



$ M ac(R,0;a) 



^(2a 2 -R 2 ) 



4o 
7rcroG 
2a 



(2a 2 - R 2 )sm-\%) + aVi? 2 - a 2 



fori? < a 
fori? > a 



(18) 



The radial force vector in the z = plane is found by taking the limit of equation [15] or by 
differentiating equation [18] with respect to R. 



FR,Mac(R,0;a) 



n 2 R<T G 

2a 
irapG 



Rsin~ l (a/R) - a^Jl - a 2 /R 2 



fori? < a 
fori? > a 



(19) 



The axial force vector on the z axis is found by taking the limit of equation [15] or differentiating 
equation [T71 with respect to z. 



F z ,Mac(0,z;a) 



2trj G 



a 



zsm 



-i 



n 



Va 2 + z 2 



a sgn(z) 



(20) 



4. The n=2 Disk 



Gonzalez fe Reinal (120061 ) give a closed form solution for the potential of the n=2 disk in elliptic 



coordinates. 

The disk surface density of the n=2 disk is 

r a (l-i? 2 /a 2 ) 3 / 2 fori?<a 



Y. D 2(R;a) 







for i? > a 



(21) 



This mass distribution can be obtained from equation [10] the disk surface density of the n=l 
disk, with the transformation: 



3 '■" 

a 3 j,, 



T> D2 (R; a) = -^j / d 2 £M ac (i?;d)<id 



(22) 



The corresponding potential is: 



3 

a 3 ./o 



$m(R, z;a) = ^ / d 2 $Mac(-R, ^;d) dd 



(23) 
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Where the expression for $Mac is given by equation [13] above. Equation [23] contains terms which 
after change of variable have with the form J x n sin~ 1 (x)dx and J x n \J\ ± x 2 dz and can be solved 
using iGradshteyn &: Ryzhikl (11994] ) 2.262, 2.813, and 2.833. After collecting terms, the result is 

_ 1 / a + I z* 



reasonably compact: 

7T(7oG 



64a 3 



3(8a 4 - 8a 2 R 2 + 16a V + 3i? 4 - 24R z z z + 8z 4 ) sin 



R 



sin 



a — I z 
~R~ 



+ a(18a 2 - 9R 2 + 26z 2 ) (V-R 2 - (a + I z) 2 + R 2 - (a - I z) 2 



(58a 2 - 55i? 2 + 50z 2 ) ( I ^/i? 2 - (a + I z) 2 - I V 'i? 2 - (a - I z) 2 



(24) 



Equation [23] can be converted to an entirely real expression by using the identities IA11 1X51 and 



<S> D2 (R,z;a) 



7R7oG 

64a 3 



6(8a 4 - 8a 2 R 2 + 16a 2 z 2 + 3i? 4 - 2AR 2 z 2 + 8z 4 )sin 



i ( f i — h 
2R 



+ V2a(l8a 2 - 9R 2 + 26z 2 ) V/1/2 - h 



(25) 



V2\z\ (58a 2 - 55i? 2 + 50z 2 ) + /s 



where /1, /2, /3 are given by equation [9] above. 

The gravitational field of the n=2 disk is the gradient of the potential using <3? as given by 
equation [25] The resulting expressions were expressed as entirely real functions by using the 
identities by using identities [AH [M] and \M\ ■ 



F R ^ D2 {R,z;a) 



37T<7oG 

16i?a 3 



2R 2 (4a 2 R 2 - 3R 4 + 12i? 2 z 2 )sin" 



2R 

, vhh^h 



V2a{2a 4 - 5a 2 R 2 + 4a 2 z 2 + 3i? 4 - 25i? 2 z 2 + 2z 4 )- 



V2\z\ (2a 4 + 9a 2 i? 2 + 4a 2 z 2 - 13i? 4 - lli? 2 z 2 + 2z 4 ) 



hh 

V/1/2 + h 



hh 



7T(ToG 

4a 3 



- 6z(2a 2 - 3i? 2 + 2z 2 )sin- 1 ( ^— ^ 

V 2R 



V2az(13a 2 - 13i? 2 + 17z 2 ) 



V/1/2 ~ /3 

hh 



V2 sgn(z)((4a 4 - 8a 2 R 2 - 3a 2 z 2 + 4i? 4 - 7i? 2 z 2 - llz 4 ) 



V/1/2 + h 
hh 



where fi, h, h are given by equation [9] above. 



(26) 



(27) 
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The potential on the z axis and on the z = plane can be found by by taking the limit of 
equation | 



$D 2 (0,z;a) 



7T(JoG 

4a 3 



3 (a 4 + 2a 2 z 2 + z 4 ) sin" 1 



a 



\/a 2 + z 2 



a | z | (5a 2 + 3z 2 ) 



(28) 



$ D2 (i?,0;a) 



-3g£oG( 8 a 4 -8a 2 i? 2 + 3i? 4 ) 



37TQ-QG 

32a 



8a 2 i? 2 + 3i? 4 )sin" 1 (§) + 3a(2a 2 - i? 2 )\/i? 2 -a 2 



fori? < a 
fori? > a 
(29) 



The radial force vector in the z = plane is found by taking the limit of equation 
differentiating equation [29] with respect to i?. 



or 



F R , D2 (R,0;a) 



37r igff G (4a 2 - 3i? 2 ) for i? < a 

[i?(4a 2 - 3i? 2 )sin- 1 (a/i?) - a(2a 2 - 3i? 2 ) y 7 ! - a 2 /# 2 l fori? > a 

(30) 

The axial force vector on the z axis is found by taking the limit of equation [57] or differentiating 
equation [28l with respect to z. 



F z ,D2(0,z;a) 



7T(7oG 

2a 3 (a 2 + z 2 ) 



at 'I i 2\2 ■ -1 

o(a + z ) zsm 



a 



Va 2 + 



(31) 



+ a^ 2 (13a 2 + 17**) + asgn(z)(4a 4 - 3a 2 z 2 - \\z 



4.1. Comparison of the Maclaurin disk and the n=2 disk 

Table 1 compares important properties of the two disks. Figure [Th. compares the surface mass 
density and potential. Figure [lb compares rotational velocity in the disks. The n=2 disk is more 
centrally concentrated than the Maclaurin disk. The rotational velocity increases more quickly in 
the inner disk and begins to fall before reaching the edge of the disk. As is apparent from figure 
[lb, the derivative of the circular velocity of the Maclaurin disk is discontinuous at the edge of the 
disk whereas the n=2 disk is better behaved. 
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Fig. 1. — Compare the surface mass density and circular velocity, V c = y (—VFr(R, z), for the 
Maclaurin disk and the n=2 disk. Reduced units are used so that the masses of both disks and the 
point mass are 1.0 



Table 1 comparison of the Maclaurin and the n=2 disks 



Property 




Maclaurin disk 


n=2 Disk 




Surface density 


£(i?) = 


aoVl " R 2 /» 2 


= a (l - R 2 /a 2 ) 3 / 2 




Total mass 


M = 


2 9 

g7ro : (To 






Circular velocity 


V 2 (R,0) = 


n 2 R 2 a G 
SitR^MG 


37r 2 R 2 a G{4a 2 - 

n 16 « 3 n 

l5TrR 2 MG(4a 2 - 


3R 2 ) 
3R 2 ) 


Disk edge velocity 


V 2 (a,0) = 


4a 3 
7T cua^G 

2 


32a 3 

3ir qo"oG 

16 
15^MG 








3vrMG 








An 


39n 





5. Example: The force field of a simple galaxy model 

Three-dimensional problems such as that of the structure and kinematics of the extra-planar 
gas will benefit from the use of the new density-potential pairs. A simple galaxy model was 
constructed for illustration. The model consists of a an n=2 disk and a core/bulge region modeled 
as a point mass. This model is defined by three parameters: the mass of the disk, the mass of 
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Fig. 2. — Circular velocity, V c = \J {—VFr(R, ~zj, of a simple galaxy model at z=0, z=2.6 kpc, 
and z=3.9 kpc above the disk. The model consists of an n=2 disk and a point mass to represent 
the core + bulge. The masses of the two components are Mdisk = 3.5 x 1O 1O M0 and Mp i n tMass = 
4 x 1O 1O M . The radius of the disk is a = 19 kpc. 



the core/bulge region, and the diameter of the disk. As shown in figure [21 the circular velocity, 
calculated as V c = \J (—VFr(R, z), is nearly constant over much of the disk. Also, the derivative 
of the circular velocity with z is nearly linear over a wide range of both R and z. 



Figure [2] agrees surprisingly well with figure 5 of iFraternali et al.l (|2005l ) which shows th at the 
measu r ed velocity of HI for NGC891 decreases li nearly with the height above t h e disk. See alsolRand 
1997 ): S waters et all d 19971) : Kamphuis et al.1 (120071 ): bosterloo et all (120071 ) : IFraternali fc Binnev 



(2006); 



Barnabe et al. ( 



20061 ). Further work is planned on this topic. 



6. Summary and conclusion 

We have presented new solutions for a family of finite disks. Closed form expressions in 
cylindrical coordinates using elementary functions are given for the potential and gravitational 
force for the disks with surface density T, n = o"o(l — R 2 / a 2 ) n ~ l l 2 with n = 0, 1, 2. Expressions are 
also given for the limiting cases of R = and z = 0. 

These solutions fill a need and should make it easier to model 3-D gravitational phenomenon 
involving disk galaxies. This is particularly important due to the recent availability of detailed 
kinematic data above the plane of the disk. 
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I am grateful to the anonymous referee for useful suggestions and comments which improved 
the presentation. 
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A. Some simple identities 

A number of identities are gathered here. In all cases x,y S 3? The ranges of validity must 
avoid the discontinuity of the principal value of the square root function on the negative real axis. 

Equation IA1I can be proved by taking the sin of both sides; reducing the terms using the 
identities sin(a + b) = sin(a)cos(6) + cos(a) sin(6) and sin(2a) = 2 sin(a) cos(6); and substituting 
cos = yl — sin 2 . The other identities can be found by substituting into the relation 



y/x + Iy 



\J x 2 + y 2 + x + sgn(y)l J yj x 2 + y 2 



V2 

and into the expression found by taking the reciprocal of both sides. 

2sin _1 



sin 1 (x — I y) - 


+- sin" 


^(x + Iy) = 


\/x - 


-iy 


+ y/x + Iy = 


ly/x- 


iy- 


ly/x + Iy = 


1 




1 




iy 


y/x-ly 


I 




I 




iy 




\/l-(x + Iyf + 


^- 


-(x-lyf = 



\^J{x + l) 2 + y 2 - y(x - l) 2 + y 2 (Al) 



\/2 y/ Va: 2 +j/2+a; (A2) 

sgn(y)v^\/ V^ 2 + y 2 - a: (A3) 

= v / 2 V - (A4) 

V £ + y 1 

sgn(y)V2 VV (A5) 

v # + y 



V2\J a/I - 2x 2 + 2y 2 + x 4 + 2a; 2 y 2 + y 4 + 1 - re 2 + y 2 (A6) 

y^l - (x + 1 y) 2 - I \Jl- {x- lyf = sgn(xy)v / 2 ^/ a/1 - 2x 2 + 2y 2 + a; 4 + 2a; 2 ?/ 2 + y 4 ~ 1 + x 2 - y 2 

(A7) 



_ . . Jl - 2x 2 + 2y 2 + x 4 + 2a; 2 y 2 + y 4 + 1 - .t 2 + y 2 
\/2— (A8) 



l-Or-Iy) 2 \/l-(x + Iy) 2 V^W+^+W+? 



_ V 1 - 2a;2 + V + xi + 2x2 y 2 + y i -i + x 2 

sgn(xy) V2 ====== ===== =— 



l-(x-lyf ^-(x + Iy) 2 V 'l-2x 2 + 2y 2 +x 4 + 2x 2 y 2 +y 4 

(A9) 
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